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Haemophilus influenzae is a host adapted human pathogen known to contribute to a 
variety of acute and chronic diseases of the upper and lower respiratory tract as well 
as the middle ear. At the sites of infection as well as during growth as a commensal the 
environmental conditions encountered by H. influenzae will vary significantly, especially 
in terms of oxygen availability, however, the mechanisms by which the bacteria can 
adapt their metabolism to cope with such changes have not been studied in detail. 
Using targeted metabolomics the spectrum of metabolites produced during growth of 
H. influenzae on glucose in RPMI-based medium was found to change from acetate as 
the main product during aerobic growth to formate as the major product during anaerobic 
growth. This change in end-product is likely caused by a switch in the major route of 
pyruvate degradation. Neither lactate nor succinate or fumarate were major products of 
H. influenzae growth under any condition studied. Gene expression studies and enzyme 
activity data revealed that despite an identical genetic makeup and very similar metabolite 
production profiles, H. influenzae strain Rd appeared to favor glucose degradation via 
the pentose phosphate pathway, while strain 2019, a clinical isolate, showed higher 
expression of enzymes involved in glycolysis. Components of the respiratory chain were 
most highly expressed during microaerophilic and anaerobic growth in both strains, but 
again clear differences existed in the expression of genes associated e.g., with NADH 
oxidation, nitrate and nitrite reduction in the two strains studied. Together our results 
indicate that H. influenzae uses a specialized type of metabolism that could be termed 
"respiration assisted fermentation" where the respiratory chain likely serves to alleviate 
redox imbalances caused by incomplete glucose oxidation, and at the same time provides 
a means of converting a variety of compounds including nitrite and nitrate that arise as 
part of the host defence mechanisms. 
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INTRODUCTION 

Haemophilus influenzae is a host-adapted human pathogen that 
occurs naturally in the nasopharynx of healthy children and 
adults (Kuklinska and Kilian, 1984). However, it is also involved 
in many acute and chronic infections such as otitis media, asthma 
and chronic obstructive pulmonary disease (COPD) (Karasic 
et al., 1989; Berkovitch et al., 2002; Revai et al, 2008; Essilfie 
et al., 2011, 2012; Moghaddam et al., 2011). The occurrence of 
these infections is linked to the ability of H. influenzae to per- 
sist in different niches in the human body, including the lung, 
nasopharynx, and the middle ear. These niches vary in car- 
bon metabolite and oxygen availability and also in ambient pH. 
For example, while the lung is a mostly aerobic environment 
(Worlitzsch et al., 2002) the middle ear is almost completely 
anaerobic (Luntz et al, 1995; Sade et al., 1995). It follows that 
H. influenzae, like other bacterial pathogens, should possess a 
repertoire of metabolic pathways that enables it to survive in the 
variety of the environments that it encounters in the human host. 

Recent studies on pathogens such as Mycobacterium species 
(Rhee et al., 2011), Legionella pneumophila (Edwards et al., 



2009, 2013), and Escherichia coli (Alteri and Mobley, 2012) have 
demonstrated clearly that subtle changes in pathway composi- 
tion or enzyme activities affect virulence of these bacteria. For 
example, Mycobacterium bovis showed attenuation in mice when 
it lacked phosphenolpyruvate (PEP) carboxykinase, an enzyme 
involved in the synthesis of PEP from oxaloacetate (Liu et al., 
2003). Similarly, the glyoxylate cycle and the TCA cycle enzymes 
isocitrate lyase and succinate dehydrogenase have been shown to 
be critical for virulence and persistence of Mycobacterium tuber- 
culosis and E. coli (McKinney et al., 2000; Alteri and Mobley, 
2012). The metabolic processes underpinning the ability of 
H. influenzae to persist in the human body have never been stud- 
ied in depth and analysis of metabolism in this bacterium has 
been limited mostly to genome reconstructions and the genera- 
tion of in silico models of metabolic fluxes (Edwards and Palsson, 
1999; Schilling and Palsson, 2000; Schilling et al, 2000; Papin 
et al, 2002; Raghunathan et al, 2004). 

Early work indicated that at least 67% of H. influenzae pro- 
teins had homologs in E. coli and that the smaller genome 
size of H. influenzae was mostly caused by loss of paralogs 
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(Tatusov et al., 1996). Reconstruction of the potential metabolic 
pathways present in the H. influenzae Rd genome (Edwards 
and Palsson, 1999) showed that this strain possesses a complete 
Embden-Meyerhof-Parnass (EMP) pathway (glycolysis) and pen- 
tose phosphate pathway (PPP) for glucose degradation, but lacks 
several enzymes of the TCA cycle (Figure 1). This latter prop- 
erty is typically associated with fermentative bacteria in which the 
TCA cycle is mostly used to produce precursors for biosynthesis. 
In H. influenzae the functional minimization of the TCA cycle is 
extreme with most enzymes of the oxidative branch being absent, 
including citrate synthase that mediates the initial condensation 
of acetyl-CoA and oxaloacetate (Figure 1). 

Several other studies have used genomic and/or flux analysis 
approaches to date, where in most cases the only allowed products 
of H. influenzae metabolism were acetate, carbon dioxide and 



in some cases certain amino acids (Edwards and Palsson, 1999; 
Schilling and Palsson, 2000; Schilling et al, 2000; Papin et al, 
2002). All of these studies appear to have aimed at determining 
the robustness of the H. influenzae metabolic network as well 
as the identification of stable flux states and the associated phe- 
notypes. A combination of proteomic data and bioinformatics 
was used by Raghunathan et al. (2004) to study anaerobic and 
microaerophilic growth of H. influenzae, and led to the identi- 
fication of formate and acetate as metabolic end-products using 
specific enzymatic detection protocols. Interestingly, the amounts 
of acetate and formate produced did not vary in response to 
changing oxygen availability (Raghunathan et al., 2004). Given 
the importance of central carbon metabolism in bacterial sur- 
vival and colonization in the host we initiated a study of carbon 
metabolism in two strains of H. influenzae, a clinical isolate from 
a COPD patient (Hi2019) (Campagnari et al, 1987) and the ref- 
erence strain Rd KW20 (HiRd) (Fleischmann et al., 1995) using 
a combination of targeted NMR-based metabolomics and qRT- 
PCR with the aim of understanding how this pathogen adapts its 
metabolism to varying oxygen availability. 

METHODS 

BACTERIAL STRAINS AND GROWTH CONDITIONS 

H. influenzae Rd KW20 (Fleischmann et al, 1995) and 2019 
(Campagnari et al, 1987) were grown at 37°C for 16-24 h using 
either brain heart infusion (BHI) broth (Difco) supplemented 
with hemin (10 (jig/ml) and NAD (2 u,g/ml) or chemically defined 
medium (CDM). 235 ml of CDM contained 191ml RPMI 1640 
(catalog no.: 11879-020; Life Technologies), 5.8 ml of 1M HEPES 
(pH 7.2-7.5, cat no.: 15630-080, Life Technologies), 2 ml of 
100 mM MEM sodium pyruvate (catalog no.: 113600-070; Life 
Technologies), 10 ml uracil (2 mg/ml), 20 ml inosine (20 mg/ml), 
2 ml NAD (1 mg/ml), 4 ml hemin (1 mg/ml), and 2.35 ml 1M 
glucose (final cone: 10 mM) as described in Coleman et al. 
(2003). BHI agar plates were supplemented with Levinthal's base 
(Alexander, 1965; Cooper et al., 2003). Cultures were grown 
under aerobic, microaerophilic, or anaerobic conditions as in 
(Cooper etal., 2003). 

GROWTH CURVES 

Bacterial growth in liquid culture was determined in triplicate 
using CDM medium and three different oxygen concentrations 
(Cooper et al., 2003). Growth experiments used sterile 250 ml 
flasks (aerobic growth: 25 ml medium, shaking at 200 rpm, oxy- 
gen transfer coefficient (A^a) 87.4 h _1 ; microaerophilic growth: 
150 ml medium, shaking at 200 rpm, oxygen transfer coefficient 
{kid) 11.5h _1 ) or completely filled 50 ml tubes. Main cultures 
were inoculated to an initial ODgoonm of 0.15 (aerobic growth) 
or 0.05 (microaerophilic and anaerobic growth) from overnight 
cultures grown from glycerol stocks. Samples for ODgoo nm deter- 
mination (0.1ml) were taken every half hour for the first 4h 
(aerobic growth) and every hour for microaerophilic and anaer- 
obic growth for the first 8 h with a final reading for all cultures 
being taken after 24 h. All growth curves used triplicate samples 
for each growth condition tested on the day they were carried out, 
and repeat experiments (biological replicates) were conducted at 
least once for each growth condition and bacterial strain. 
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FIGURE 1 | Schematic representation of H. influenzae pathways for 
central carbon metabolism (A) and the elements of the respiratory 
chain (B). Abbreviations: pyruvate dehydrogenase complex, aceF; Pyruvate 
formate lyase, pfIA; Acetate kinase, ackA; NAD+-dependent D-lactate 
dehydrogenase, IdhA; Formate dehydrogenase, fdxG; NADH 
dehydrogenase, ndh; NADH dehydrogenase, nqr; L- Lactate 
dehydrogenase, lldD; D-Lactate dehydrogenase, did; Cytochrome bd 
oxidase, cydA, cydB; DMSO reductase, dmsA; Nitrate reductase, napA; 
TMAO reductase, torZ; Nitrite reductase, nrfA; Fumarate reductase, frdA; 
malate dehydrogenase, mdh; PEP carboxylase, ppc; PEP carboxykinase, 
pckA; 2-oxoglutarate dehydrogenase, sucAB; Succinyl-CoA synthetase, 
sucCD; fumarate hydratase, fumC. 
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PREPARATION OF NMR SAMPLES AND NMR SPECTROSCOPY 

Samples for ID H NMR were collected at three different 
time points, corresponding to early, mid and late exponential 
growth phases. Collection time points were selected based on 
growth curve experiments conducted in preparation for this 
work. To quantify metabolite concentrations in the medium 
0.9 ml samples of H. influenzae culture were pelleted by centrifu- 
gation (3min, RT, rcf: 20,000), the supernatants were removed 
to a fresh tube followed by the addition of azide to a final 
concentration of 0.01% from a 2% (w/v) stock and storage 
at — 80°C. Samples for *H NMR analysis had a volume of 550 (xL 
and contained 500 u,L of preserved culture supernatant, 10% 
D2O, 18 (xM 2,2-dimethylsilapentane-5-sulfonic acid (DSS) and 
18 (jlM 1,1-difluoro-l-trimethylsilanylphosphonic acid (DFTMP) 
in 5 mm NMR tubes. Proton NMR spectroscopy was performed 
on a Bruker AV500 500 MHz spectrometer (Bruker Biospin, 
Rheinstetten, Germany), with a 5 mm self-shielded z-gradient 
triple resonance probe. One-dimensional NOESY spectra were 
measured as previously described (Schirra et al, 2008) at 298 K 
with 128 scans at 32 k resolution with a spectral width of 14 
ppm, using the noesylprd pulse program (Bruker pulse program 
library). The water signal was suppressed with low-power contin- 
uous wave irradiation during the NOESY mixing time of 100 ms 
and the relaxation delay of 2.3 s. ID spectra were processed using 
Topspin 3.0 (Bruker Biospin). Spectra were manually phased and 
baseline corrected, and chemical shifts were referenced to the DSS 
signal. Metabolites were identified by comparing putative identi- 
fications obtained using the Chenomx NMR Suite 7.1 (Chenomx 
Inc., Edmonton, Canada) with 'H chemical shifts from either 
1 mM pure metabolite standards run separately, or by comparison 
with data from public databases such as BioMagResBank [www. 
bmrb.wisc.edu] or the Human Metabolome Database [www. 
hmdb.ca/]. NMR signals of metabolites of interest were man- 
ually integrated using Topspin 3.0 (Bruker Biospin) to obtain 
absolute integral intensities. Absolute integrals were calibrated in 
Chenomx NMR Suite 7. 1 against the intensity of the DSS signal 
(18 |xM) in the NMR spectrum of the uninoculated medium to 
obtain metabolite concentrations. Concentrations of metabolites 
not present in the uninoculated medium (labeled "Products" in 
Tables S5, S6) were calibrated relative to their concentration in 
the sample of strain 2019 grown for 4h under micro aerophilic 
conditions. NMR data of the metabolites of interest and the 
NMR signals used for concentration determination are listed in 
Table SI. 

QUANTITATIVE RT-PCR 

Standard techniques were used throughout (Ausubel et al, 2005). 
RNA was isolated from 1.5 ml to 5 ml of H. influenzae cul- 
tures grown to an ODgoo of 0.4. Culture samples were either 
pelleted by centrifugation followed by immediate freezing at 
— 80°C or preserved using RNAprotect Bacteria Reagent (Qiagen) 
prior to freezing at — 80°C. RNA was extracted using either the 
RNeasy Mini Kit (Qiagen) or the GE Healthcare Illustra RNAspin 
Mini Kit according to the manufacturer's instructions. RNA 
samples were tested for gDNA contamination using standard 
PCR (GoTaqgreen Mastermix, Promega) with RT_16S primers 
(Table S2). If no PCR product was observed after 34 cycles, 



samples were considered to be gDNA-free. RNA concentrations 
were determined using the Quant- it RNA kit (Life Technologies). 
cDNA was synthesized from 500 ng of RNA using the Superscript 
III Reverse Transcriptase (Life Technologies) and RNAsin RNase 
inhibitor (Promega). qRT-PCR (10 |xl) reactions used diluted 
cDNA (1:10-1:100 000) as template, the SYBR green Mastermix 
(Applied Biosystems) and primers (Table S2) designed in Vector 
NTI Advance (Life Technologies) to produce PCR products of 
120 bp. The 16S rDNA gene was used as the reference gene. For 
each individual experiment, qRT-PCR reactions were set up in 
triplicate (technical replicates) using an epMotion workstation 
(Eppendorf) and run on an ABI 7900 sequence detector. qRT- 
PCR experiments targetting all monitored genes were repeated at 
least three times for each strain with similar result being obtained. 
Data analysis and normalization was performed as in (Kappler 
etal, 2005). 

ENZYME ASSAYS 

Cell cultures grown in BHI supplemented with 20 mM bicar- 
bonate for 16-24 h at three different oxygen tensions were har- 
vested by centrifugation (2370 x g, 4°C, 15min). Cell pellets 
were stored at — 20° C until further use. Cell extracts were pre- 
pared by resuspending the frozen pellets in 2 ml BugBuster 
Mastermix (Novagen). The samples were then incubated on a 
shaker at room temperature for 20 min followed by centrifuga- 
tion (20238 x g, 5 min, RT). The supernatants were removed 
to clean tubes and were used for enzyme assays. DMSO reduc- 
tase was assayed at 37°C as in (Jones and Garland, 1977) 
using a Hitachi UV3000-Spectrophotometer. The reaction was 
monitored at 600 nm using the oxidation of benzyl viologen 
(e 60 o = 7.4 mM -1 cm -1 ) in the presence of either 17 mM DMSO 
or 5 mM MetSO (Lester and Demoss, 1971). Formate dehydroge- 
nase activity was assayed according to (Enoch and Lester, 1982) by 
monitoring the reduction of dichlorophenolindophenol (DCPIP; 
66oo = 21mM~ 1 cm -1 ) (Lester and Demoss, 1971). All assay 
mixtures were degassed with nitrogen for 2 min before addition 
of the reductant and the cell extracts and used cuvettes sealed 
with rubber septa (Sigma Aldrich). Assays were conducted using 
at least two biological replicates (grown on different days) and at 
least three technical replicates (repeat assays) for each cell extract. 
Protein concentrations were determined using the Sigma- Aldrich 
BCA- 1 protein determination kit. 

RESULTS AND DISCUSSION 

CENTRAL CARBON METABOLISM IN H. influenzae — PATHWAYS AND 
ENERGY GENERATION 

Analysis of available H. influenzae genomes clearly showed 
that the organization of the central carbon metabolic pathways 
(Figure 1) (Table S3) is highly conserved and confirmed results 
from previous work that had focused solely on the H. influenzae 
reference strain Rd (Fleischmann et al., 1995). All H. influenzae 
genomes encode enzymes for the EMP pathway (glycolysis), with 
the exception of glucokinase, and contain all enzymes of the PPP 
Pyruvate, the main product of hexose degradation, can then be 
converted into a variety of other compounds (Figure 1). Acetyl - 
CoA can be produced using pyruvate dehydrogenase (encoded 
by aceEF, IpdA) or pyruvate formate lyase (pflA). Pyruvate 
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can also be reduced to D-lactate by NAD-dependent lactate 
dehydrogenase (LdhA). 

Of the TCA cycle enzymes only malate dehydrogenase, 
fumarate hydratase, succinyl CoA synthetase, and ot-ketoglutarate 
dehydrogenase are present. All of these enzymes are important for 
the generation of biosynthetic precursors and the remnant TCA 
cycle is also linked to the respiratory chain via fumarate reduc- 
tase. An exception to this is H. influenzae strain PittGG where 
the frdA gene encoding the large subunit of fumarate reduc- 
tase contains a frameshift mutation, suggesting loss of function. 
Carbon units can enter this partial TCA cycle via the reaction of 
phosphenolpyruvate carboxykinase while acetyl-CoA produced 
by pyruvate dehydrogenase or pyruvate formate lyase can be 
converted to acetate via acetyl-P producing one ATP (Figure 1). 

Incomplete TCA cycles are typically found in fermentative 
bacteria but in contrast to these, H. influenzae strains possess a 
versatile respiratory chain with five dehydrogenases (including 
NADH, lactate and formate dehydrogenases) transferring elec- 
trons into the menaquinone (MQ) pool, and at least six terminal 
reductases transferring the electrons to a variety of electron accep- 
tors including oxygen, dimethylsulfoxide (DMSO) and various 
nitrogen compounds (Figure 1, Table S3). The components of 
this respiratory chain are almost completely conserved [excep- 
tion: formate dehydrogenase and lactate dehydrogenase (LldD) 
in strain PittAA; torZ in strains 10810, 86-028NP, F3031, and 
22.4-21]. This allows H. influenzae strains to use this versatile res- 
piratory chain to easily re-oxidize the NADH produced by the 
partial oxidation of carbon sources under a variety of growth 
conditions. 

We extended our analyses to include available genomes of 
other Haemophilus species (Table S4) which showed that virtu- 
ally all Haemophilus genomes encode both glycolysis and PPP 
pathways. The presence of a partial TCA cycle is limited to H. 
influenzae, Haemophilus parainfluenzae, Haemophilus haemolyti- 
cus, and Haemophilus ducreyi, while most Haemophilus parasuis 
and Haemophilus somnus species appear to contain a complete 
TCA cycle. All Haemophilus species appeared to have reason- 
ably versatile respiratory chains although some features such 
as the presence of the dmsA and torZ terminal reductases as 
well as three lactate dehydrogenases and especially the UdD and 
did encoded enzymes, appear to be specific to H. influenzae 
(Tables S3, S4). 

METABOLIC END PRODUCTS PRODUCED BY H. influenzae STRAIN RD 
(HIRD) CHANGE IN RESPONSE TO CHANGING OXYGEN TENSIONS 

In keeping with the observation that H. influenzae only possess 
an incomplete TCA cycle, organic acids such as acetate and for- 
mate were the main products detected by ^-NMR following 
aerobic, microaerophilic and anaerobic growth of H. influenzae 
in chemically defined medium (CDM) with glucose (Table 1). 
The two main carbon sources in the medium were differentially 
utilized. While pyruvate was quickly consumed (within 4h fol- 
lowing inoculation) glucose was initially produced by aerobic and 
microaerophilic cultures, presumably by gluconeogenesis, before 
being utilized in mid to late exponential phase (Figure 2). Acetate 
was the major end product during aerobic growth (Tables 1, S5), 
indicating that the dominant route of pyruvate catabolism was 



generation of acetyl-CoA and further conversion to acetate via 
acetyl phosphate. 

Under microaerobic conditions formate started to accumulate 
in addition to acetate (Tables 1, S5, Figure 2). This likely repre- 
sents a switch in pyruvate metabolism from pyruvate dehydroge- 
nase to an increased flux through pyruvate formate lyase (PFL), 
which during anaerobic growth of H. influenzae would lead to 
the observed accumulation of formate as the major end-product 
(Tables 1, S5). 

Despite being typical products of bacterial "mixed acid" type 
fermentations, small organic acids such lactate and succinate 
were only produced in small amounts under all three conditions 
studied here (Tables 1, S5). Other metabolic endproducts were 
hypoxanthine (0.7-3.5 mM), probably resulting from the degra- 
dation of inosine, glycerol and glycine (4.38-1.04 mM, amounts 
decreasing with decreasing oxygen tensions) (Tables 1, S5). 

Our results confirm that acetate and formate are the main 
metabolic end products of H. influenzae metabolism, but in con- 
trast to a previous study that reported little or no variation 
in acetate (0.46-0.6 mM) and formate (0.29-0.3 mM) concen- 
trations in microaerophilic and anaerobic H. influenzae cul- 
tures (Raghunathan et al, 2004), our data demonstrate a clear 
shift from acetate to formate as the main metabolic product 
in response to decreasing oxygen availability. The difference in 
the observed endproduct concentrations may be due to the use 
of different growth media. While this study used a chemically 
defined medium with glucose as the main carbon source, the 
study by Raghunathan et al. (2004) used a complex medium, 
Brain-Heart Infusion, that would likely have had a higher propor- 
tion of peptides and amino acids and a lower sugar contents than 
the tissue culture medium based CDM used in our work. While 
Brain-Heart Infusion is a standard medium for the cultivation 
of H. influenzae, the CDM formulation should more accurately 
reflect the conditions encountered by the bacteria when grown 
in contact with human cells. Preliminary experiments also (data 
not shown) indicated that the matrix of Brain-Heart Infusion is 
unsuitable for NMR analyses. 

Our observations of changes in the relative amounts of 
metabolic endproducts produced by H. influenzae are consis- 
tent with constraints of the ability of the bacteria to maintain 
redox homeostasis in the absence of alternative external electron 
acceptors once oxygen becomes limiting. During aerobic growth 
when acetate is the main metabolic product NADH generated 
during glucose breakdown can be reoxidized by the respiratory 
chain and e.g., the cytochrome bd oxidase. Once oxygen becomes 
limiting, there is a switch to production of formate using PFL, 
which does not generate NADH. The alternative would be for 
pyruvate dehydrogenase to continue to be the major enzyme 
of pyruvate catabolism with NADH being consumed via lac- 
tate dehydrogenases and fumarate respiration, but the low lev- 
els of D-lactate and succinate produced suggest that these are 
minor pathways. 

METABOLIC ADAPTATION IS SUPPORTED BY UNDERLYING CHANGES 
IN GENE EXPRESSION 

19 genes encoding key enzymes involved in H. influenzae carbon 
metabolism and respiration were selected (Table S2) and changes 
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Table 1 | Key metabolites detected in the growth medium of H. influenzae strain RD. 



Metabolite 


Initial cone (mM) 


Aerobic (mM) 


Microaerophilic (mM) 


Anaerobic (mM) 


Substrates Glucose 


10 


8.326 


4.228 


4.082 


Pyruvate 


0.87 


0.031 


0.040 


0.025 


Inosine 


6.5 


0.866 


2.444 


2.689 


Products Formate 


n.d. + 


1.86 


4.05 


23.36 


Acetate 


n.d. + 


11.61 


6.43 


4.28 


Hypoxanthine 


n.d.+ 


3.53 


1.75 


0.701 


Glycerol 


n.d.+ 


0.856 


0.583 


0.185 


Succinate 


n.d.+ 


0.011 


0.034 


0.641 


Lactate 


n.d.+ 


0.003 


0.036 


0.013 



Samples were taken in late exponential phase. A list of minor metabolites detected is available in Table SS. + not detected. 








2019 glucose 











FIGURE 2 | Changes in metabolite concentrations in cultures of H. 
influenzae RD and 2019 grown under aerobic (closed circles), 
microaerophilic (open circles), and anaerobic (closed triangles) 
conditions. The data shown is for two growth substrates, pyruvate and 
glucose (rows 1 and 2) and two key metabolites produced during growth, 
acetate, and formate (rows 3 and 4). 



in gene expression determined under the same conditions used 
for the NMR experiments. 

Interestingly, under all growth conditions tested high levels 
of pyruvate dehydrogenase (aceF), acetate kinase (ackA), and 
the respiratory L-lactate dehydrogenase (UdD) were detected, 



indicating key roles for these enzymes in H. influenzae 
metabolism (Figure 3). In fact, HiRd is capable of growth using 
lactate as a carbon source (data not shown), and the LldD enzyme 
could be important under those conditions. It has also been sug- 
gested that LldD might be important during invasive disease as 
lactate is commonly present in blood (Wong et al, 2007). 

Based on the metabolic model (Figure 1) the initial degra- 
dation of glucose in H. influenzae would occur using either 
the PPP or glycolysis pathways, and interestingly both 6- 
phosphofructokinase (pfkA) and glucose-6-phosphate 1- dehy- 
drogenase (zwf) which represent glycolysis and PPP, respectively, 
as well as pyruvate kinase (pykA) were expressed most highly 
under microaerophilic and anaerobic conditions. However, rel- 
ative expression levels of zwf (and pykA) were much higher than 
those ofpfkA (16-27 fold and 24-31 fold, respectively) (Figure 3) 
suggesting that HiRd uses mainly the PPP for glucose metabolism. 
Two enzymes that would use the product of glucose break- 
down, pyruvate, pyruvate-formate lyase (pflA), and the NAD + - 
dependent LdhA lactate dehydrogenase both showed increased 
expression levels during microaerophilic and anaerobic growth. 

Of the respiratory chain enzymes all five respiratory dehy- 
drogenases were most highly expressed under microaerophilic 
conditions, with the two lactate dehydrogenases, LldD and 
Did and the sodium-translocating NADH dehydrogenase (Nqr) 
being expressed at high levels, while formate dehydrogenase and 
the non-energy conserving NdhA NADH dehydrogenase were 
expressed at low levels at all times (Figure 3). Expression of nqr 
genes was also high under anaerobic conditions. 

In contrast, for all six terminal reductases including fumarate 
reductase (frdA) and the cytochrome bd oxidase relative expres- 
sion levels were highest under anaerobic conditions. The expres- 
sion pattern of the cytochrome bd oxidase is thus similar to what 
has been seen in E. coli (Cotter et al, 1990). The gene encoding 
the NapA nitrate reductase was the only one that showed similar 
expression levels under microaerophilic and anaerobic condi- 
tions, while genes encoding sulfoxide reductase related enzymes 
(dmsA, torZ) and nitrite reductase (nrf) were only expressed at 
significant levels under anaerobic conditions. This is consistent 
with the role of nitrate and sulfoxide reductases as alternative ter- 
minal reductases that has been described for a variety of bacteria 
(Richardson et al., 2001; McCrindle et al, 2005). 
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FIGURE 3 | Expression of genes involved in central carbon metabolism 
and respiration in H. influenzae RDKW20 for cultures grown under 
aerobic (A) microaerophilic (B) and anaerobic (C) conditions. Black bars: 
Genes relevant for glucose catabolism. White bars: Genes relevant for 
pyruvate conversions, gray bars: genes relevant for respiratory metabolism. 
Abbreviations of gene names: Glucose 6-phosphate dehydrogenase, zwf; 
pyruvate dehydrogenase complex, aceF; Pyruvate formate lyase, pfIA; 
Acetate kinase, ackA; NAD + -dependent D-lactate dehydrogenase, IdhA; 
Formate dehydrogenase, fdxG; NADH dehydrogenase, ndh; NADH 
dehydrogenase, nqrB; L-Lactate dehydrogenase, lldD; D-Lactate 
dehydrogenase, did; Cytochrome bd oxidase, cydA, cydB; DMSO 
reductase, dmsA; Nitrate reductase, napA; TMAO reductase, torZ; Nitrite 
reductase, nrfA; Fumarate reductase, frdA. 



A number of the genes studied here have been previously 
identified as being under the control of major regulators of aer- 
obic/anaerobic transitions in bacterial metabolism namely FNR 
(relevant target genes: dmsA, napA, nrfA) and the ArcAB two 
component system (relevant target genes: fdx, frd, lldD, ndh) (De 
Souza-Hart et al., 2003; Wong et al, 2007; Harrington et al, 
2009). While the three FNR regulated genes studied here showed 
the expected increase in gene expression under anaerobic condi- 
tions already reported previously (Harrington et al., 2009), the 
case for the ArcA regulated genes is not so clear cut. A study using 
a proteomic approach observed and Arc regulon of ~40 proteins 
in H. influenzae, and identified three of these using mass spec- 
trometry, namely the Fdx formate dehydrogenase and the Frd 
fumarate reductase, both of which were induced under anaerobic 
conditions, and the UdD lactate dehydrogenase that was sub- 
ject to Arc-based repression under both aerobic and anaerobic 
growth conditions (De Souza-Hart et al., 2003). These observa- 
tions largely agree with our data reported above, where fdx and 
frd are expressed at higher levels when oxygen becomes limiting, 
although for fdx we observed peak expression and activity under 
microaerophilic rather than fully anaerobic conditions. A study 
monitoring gene expression changes in an H. influenzae arcA 
mutant under anaerobic conditions identified a total of 23 genes 
controlled by ArcA which included fdx, UdD, and ndh (Wong 
et al., 2007). All three of these genes were reported to be negatively 
regulated by ArcA under anaerobic conditions which is in agree- 
ment with our gene expression data (Figure 3). The frd fumarate 
reductase was not identified as a target of Arc-regulation by 
(Wong et al., 2007). Data from our work and these previous stud- 
ies may indicate that regulation of some enzyme such as the Fdx 
formate dehydrogenase is complex and may involve more than 
one regulator. 

In summary, our data presented above indicate that H. influen- 
zae metabolism is likely adapted to conditions where low levels of 
oxygen are present as most of the available respiratory pathways, 
including oxygen based respiration, are highly expressed under 
these conditions. In H. influenzae RD there also appears to be a 
preference for glucose degradation via the PPP rather than via 
glycolysis, which allows an increased production of NADPH and 
pentose sugars for the synthesis of nucleotides. 

As expression of the formate dehydrogenase encoding genes 
appeared to be extremely low (Figure 3) we conducted enzyme 
assays for formate dehydrogenase and, as a control, DMSO reduc- 
tase to confirm the qRT-PCR results (Figure 4). In keeping with 
the qPCR results the highest level of formate dehydrogenase activ- 
ity was detected in cell grown under microaerophilic conditions 
(~350mU/mg), while DMSO reductase activity was highest in 
extracts made from anaerobically grown cells (~1500 mU/mg), 
which is in keeping with the known FNR-based reglation of 
the dmsA promoter (Harrington et al., 2009). In E. coli, both 
enzymes are known to be active under conditions when oxygen 
concentrations are low (Gladyshev et al, 1996; McNicholas et al., 
1997). 

The high expression levels of formate dehydrogenase under 
microaerophilic conditions and the high levels of formate 
detected as a metabolic endproduct during anaerobic growth 
suggest a connection between the relative activities of pyruvate 
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FIGURE 4 | Enzyme activities of DMSO reductase (A) and Formate 
dehydrogenase (B) in H. influenzae RD under aerobic, microaerophilic, 
and anaerobic growth conditions. Assays were conducted with at least 
three repetitions, the experimental error is reported as standard error of the 
mean. 



formate lyase (pflA) and formate dehydrogenase (fdxG) which 
produce and degrade formate, respectively. The presence of for- 
mate dehydrogenase activity under microaerophilic growth may 
lead to an underestimation of the actual flux of carbon via PFL 
(Figure 2, Tables 1, S5) as the formate produced would be fur- 
ther degraded to carbon dioxide. Formate dehydrogenase is an 
energy-conserving Mo enzyme in the respiratory chain and in 
combination with the cytochrome bd oxidase can generate a 
proton motive force. The results infer that under microaerobic 
conditions respiration with formate as an electron donor and oxy- 
gen as electron acceptor allows additional energy conservation via 
an increased proton motive force. Under anaerobic conditions 
and in the absence of other external respiratory electron accep- 
tors formate cannot be further oxidized since H. influenzae does 
not possess formate-hydrogen lyase. 

DO ALL H. influenzae STRAINS SHOW THE SAME PATTERN OF 
METABOLIC ADAPTATION? 

While being widely used in research, the Rd reference strain has 
been in laboratory culture for decades, and we therefore extended 
our analyses to a non-typeable strain of H. influenzae, strain 2019 
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FIGURE 5 | Expression of genes involved in central carbon metabolism 
and respiration in H. influenzae 2019 for cultures grown under aerobic 
(A) microaerophilic (B) and anaerobic (C) conditions. Black bars: Genes 
relevant for glucose catabolism, White bars: Genes relevant for pyruvate 
conversions, gray bars: genes relevant for respiratory metabolism. 
Abbreviations of gene names: Glucose 6-phosphate dehydrogenase, zwf; 
pyruvate dehydrogenase complex, aceF; Pyruvate formate lyase, pflA; 
Acetate kinase, ackA; NAD+-dependent D-lactate dehydrogenase, IdhA; 
Formate dehydrogenase, fdxG; NADH dehydrogenase, ndh; NADH 
dehydrogenase, nqrB; L-Lactate dehydrogenase, lldD; D-Lactate 
dehydrogenase, did; Cytochrome bd oxidase, cydA, cydB; DMSO 
reductase, dmsA; Nitrate reductase, napA; TMAO reductase, forZ; Nitrite 
reductase, nrfA; Fumarate reductase*, frdA. 
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(Hi2019) which is an isolate from a patient suffering from COPD 
(Campagnari et al., 1987). 

Overall, the spectrum of metabolic end products was the same 
in both H. influenzae strains, but there were subtle changes in the 
amounts as well as types of minor products present (Tables 2, S6), 
and in the consumption of substrates. For example, in Hi2019 
depletion of pyruvate was faster that in HiRD, and the final 
concentrations of glucose present in the growth medium were 
lower than in HiRD indicating a more complete utilization of 
the carbon source (Figure 2). The key metabolites produced by 
Hi2019 were similar to those observed previously, although larger 
amounts of formate (13.16 vs. 4.05 mM for HiRD) were pro- 
duced under microaerophilic conditions (Tables 2, S6). No sig- 
nificant amounts of glycerol or glycine were present in the growth 
medium following growth of HI2019 although both metabo- 
lites were produced transiently during growth of the strain, 
with glycine concentrations reaching a maximal concentration 
of ~L2 mM (data not shown). 

We then conducted studies of gene expression in H. influen- 
zae 2019 using the same set of genes used for HiRd previously 
(Figure 5). Most of the gene expression patterns were similar 
to HiRD, some significant differences were also observed. The 
expression levels of pfkA (glycolysis) exceeded those of zwf 
(PPP) between 50 and 89 times, indicating that in Hi2019 gly- 
colysis is the main pathway for glucose catabolism. The ndhA 
gene was only expressed at low levels in HiRD but is always 
expressed at high levels in Hi2019. This is interesting as the NdhA 
protein oxidizes NADH without energy conservation via pro- 
ton translocation which would have implications for the energy 
balance. Expression of the second NADH dehydrogenase com- 
plex (nqrB) appeared to follow a similar expression pattern 
as in HiRD. 

Of the other respiratory chain components two of the termi- 
nal reductases showed clearly changed gene expression patterns 
in Hi2019, namely the nrfA (nitrite reductase) and the napA 
(nitrate reductase) genes both of which seemed to be expressed 
at much higher levels than in HiRd relative to other enzymes of 
the respiratory chain. 

These data lead us to propose that despite their nearly identi- 
cal genetic makeup glucose metabolism in HiRd and Hi2019 uses 
different main pathways. In HiRd the dominant catabolic route 



is propose to involve the PPP which generates NADPH (used in 
biosynthesis rather than being used in the respiratory chain) and 
glyceraldehyde 3-P. In contrast, Hi2019 appears to use glycoly- 
sis as the main pathway, where an additional molecule of ATP 
is consumed during activation of glucose. Hi2019 may also have 
a greater dependence on the respiratory chain, especially NADH 
dehydrogenase and nitrate reductase. 

The data presented above also clearly show that despite the 
differences in the regulation and utilization of different res- 
piratory chain complexes, both strains of H. influenzae use a 
type of metabolism that could be called "respiration-assisted" 
fermentation, i.e., a fermentative type of metabolism in which 
problems with redox homeostasis due to overreduction of the 
cellular NAD + pool which can be caused by purely fermentative 
growth can be alleviated through the activities of the respiratory 
enzymes. 

Despite the presence of a respiratory chain H. influenzae 
metabolism is configured so that simple chain carboxylic acids 
will always be the main product regardless of the presence or 
absence of oxygen. As a result the regulatory mechanisms at play 
in H. influenzae can be expected to be quite different from those 
found e.g., in E. coli and other bacteria capable of aerobic respi- 
ration in the presence of a complete TCA cycle. Thus, the switch 
from acetate as the main product to formate in H. influenzae may 
be due to a change in the main route of pyruvate metabolism, 
the underlying mechanisms are likely unrelated to the "acetate 
switch" described for E. coli (Wolfe, 2005). 

An interesting feature of the H. influenzae respiratory chain is 
that it contains a surprisingly large number of components not 
associated with energy conservation via a proton gradient (e.g., 
NdhA, Nqr, Nap, TorZ) which may also represent an adaptation to 
the host environment. Differential regulation of respiratory chain 
composition and pathways associated with energy generation 
have been shown to be associated with virulence in Escherichia 
coli infecting the urinary tract (Alteri et al, 2011) and the possi- 
bility that similar mechanisms exist in H. influenzae should be 
investigated in the future. While a metabolic arrangement that 
has low efficiency in energy conservation may appear unusual, 
there is some evidence for Salmonella enterica sv typhimurium 
that an incomplete TCA cycle enhances survival of the bacteria 
in macrophages (Bowden et al., 2010) and similar effects may be 



Table 2 | Key metabolites detected in the growth medium of H. influenzae 2019. 



Metabolite 


Initial cone. (mM) 


Aerobic (mM) 


Microaerophilic (mM) 


Anaerobic (mM) 


Substrates Glucose 


10 


1.98 


2.26 


3.864 


Pyruvate 


0.87 


0.018 


0.015 


0.041 


Inosine 


6.5 


2.60 


3.043 


3.347 


Products Formate 


n.d. + 


2.09 


13.16 


23.04 


Acetate 


n.d. + 


11.17 


8.06 


4.74 


Hypoxanthine 


n.d. + 


1.63 


1.53 


1.242 


Glycerol 


n.d. + 


0.046 


0.04 


0.032 


Succinate 


n.d + 


0.007 


0.401 


1.009 


Lactate 


n.d.+ 


0.069 


0.098 


0.456 



Samples were taken in late exponential phase. A list of minor metabolites detected is available in Table S6 + not detected. 
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relevant in H. influenzae but will require further investigation. 
Interestingly, for H. influenzae serotype b it has been reported that 
a loss of either pyruvate dehydrogenase (aceEF) or a-ketoglutarate 
dehydrogenase (sucAB) leads to a significant attenuation of viru- 
lence in a mouse model of systemic disease (Herbert et al., 2003). 
This underlines the importance of metabolism and metabolic 
enzymes for H. influenzae pathogenesis and survival in vivo. 
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